Unveiling the vortex-glass phase transition in superconducting Ba(Feo.9Coo.i)2As2 
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Measurements of magnetic ac susceptibility were performed on a single crystal of superconducting 
Ba(Feo.9Coo.i)2As2 close to the conditions of optimal doping. The extremely high quality of the 
investigated sample allows the detection of dynamical-scaling behaviour associated with a glass 
transition for the flux lines in the limit of weak degree of quenched disorder. The corresponding 
critical exponent and glass critical temperature are precisely quantified for different values of the 
polarizing magnetic field up to 9 Tesla. The field-dependence of the critical exponent rules out any 
dimensionality crossover for the vortex state in the investigated range, suggesting instead a gradual 
transition among qualitatively different types of glassy phases. In particular, a Bragg-like glass is 
likely to be gradually disordered towards a more amorphous state with increasing the magnetic field. 

PACS numbers: 74.25.Uv, 74.25.Wx, 74.70.Xa 



The possibility of exploring the physical properties of 
the mixed Shubnikov phase in type-II superconductors 
over a wide range of external thermodynamical param- 
eters has lead to the most important achievements of 
the research on high-Tc materials The physics of 

vortices (or fluxoids, each one bringing a quantized unit 
$0 = hc/2e of magnetic field flux) is indeed of great in- 
terest on the level of fundamental research, where the 
system is typically modeled as a statistical ensemble of 
fllaments interacting both among themselves and with 
the defects of the crystalline environment P, [1] ■ In this 
respect, a strongly debated topic concerns the actual oc- 
currence of a real phase transition between different re- 
gions in the magnetic field - temperature {H — T) phase 
diagram of the superconducting state where dissipation 
processes occur or not 0]. An impressive amount of ex- 
perimental data, especially from transport measurements 
of current- voltage {I—V) characteristics, have been inter- 
preted in terms of scaling arguments strongly suggesting 
that the scenario of critical processes is holding (for an 
early review about scaling in high-Tc cuprate supercon- 
ductors see Ref. 

The discovery of high-Tc type-II superconductivity in 
Fe-based pnictides ,^,,6;] has recently renewed the interest 
for the investigation of the vortex state. These materials 
are known to possibly bridge the gap between high-Tc 
cuprates and low-Tc alloys. Besides Tc values interme- 
diate with respect to what is typically reported for the 
aforementioned two classes of materials, in fact, pnic- 
tides join both high upper critical magnetic fields i?c2 
and almost isotropic superconducting properties quanti- 
fied by low values for the anisotropy parameter 7. Em- 
blematic compounds in this respect are those belonging 
to the so-called 122 family, displaying /Lto-ffc2 as high as 
50 Teslas while 7 as low as 2 Q. Several experimen- 
tal results have been reported in the literature about the 
detailed features of vortex pinning in connection with 
the flux motion, the positioning of the irreversibility line 



on the H — T phase diagram and the estimates of crit- 
ical currents (see, for instance, Refs. Neverthe- 
less, to the best of our knowledge, the data have been 
rarely examined in terms of the critical behaviour as- 
sociated with a phase transition to a glassy-like state 
for the vortices, but rather in terms of thermally acti- 
vated processes in the strong-disorder limit. Critical- 
ity is in fact expected to be hidden by high degrees of 
quenched disorder in the system, including the extrinsic 
effects of boundaries between different grains in powder 
samples. In this Letter, we report on ac susceptibility 
measurements performed on an extremely high-quality 
single crystal of Ba(Feo.9Coo.i)2As2. Results are quali- 
tatively and quantitatively different from what has been 
observed from similar measurements performed on other 
pnictide superconductors and strongly hint that the sce- 
nario of a vortex-glass phase transition is at work. The 
extremely low degree of quenched disorder in the inves- 
tigated sample, moreover, allows to infer interesting in- 
sights into the detailed evolution of the glassy-phase to- 
wards more amorphous conflgurations upon increasing 
the static magnetic field. 

The investigated single crystal is of rectangular shape 
with dimensions 2a = 1.7 ± 0.025 mm, 2b = 3.0 ± 0.025 
mm and 2c = 0.05 ± 0.025 mm (where the height 2c 
is parallel to the c crystallographic axis), following the 
notation reported in Ref. JA,. Details about the crystal- 
growth procedure and preliminar characterization can be 
found in Ref. [3 According to our resistivity data, the 
sample is in conditions of almost-optimal doping with a 
value for the critical temperature in zero magnetic field 
estimated as Tc(0) ~ 23.5 K. Measurements of magnetic 
ac susceptibility Xac as a function of T were performed 
by means of a Physical Properties Measurement System 
(PPMS, Quantum Design) allowing to apply polarizing 
static magnetic fields fioHdc as high as 9 T. Measure- 
ments were always performed by slowly warming the sam- 
ple after a field-cooled (FC) procedure across Tc. The 
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amplitude of the alternating magnetic field iJac was kept 
fixed to f Oe during all the measurements while the work- 
ing frequency v^a was swept between f and 10^ Hz. Both 
H^c and iJac were applied parallel to the crystallographic 
c-axis (within an accuracy of ± 1°). For this geometry, 
the demagnetization factor is estimated as Dm — 0.95 
(see Ref. [ij). The intrinsic susceptibility Xi is then ob- 
tained from the measured one Xm following the relation 
1/Xi = 1/Xm — 47rl?in (where volume units have been 
assumed for both Xi a-nd Xm)- From now on, only the 
demagnetization-corrected Xi data will be considered and 
the subscript i will be dropped for the sake of clarity. 

Fig. [T] shows a step-like diamagnetic response of the 
real part of the magnetic ac susceptibility xLc whose po- 
sition is strongly dependent on the value of H^c , in agree- 
ment with similar results obtained on both cuprate and 
pnictide superconductors 13) U , IB" 18 1 ■ Small values for 
the transition width A are detected independently on the 
value of TJdc, namely A(90kOe) ~ 2A(250Oe) ~ 0.4 K. 
This striking sharpness hints towards the excellent qual- 
ity of the sample, ensuring that the homogeneity of the 
Co-doping is extremely high. The sample is then an opti- 
mal ground for investigating the intrinsic features of the 
superconducting phase and of the motion of the vortex 
lines. Such drop in Xaci fact, is known to separate re- 
gions inside the thermodynamical superconducting phase 
where dissipative processes associated with the motion of 
vortex lines take place or not. In particular, due to the 
complex interplay of several factors like the mutual re- 
pulsion among vortices, the decreased amount of thermal 
energy and the interaction with lattice defects acting like 
pinning centers, in the low-T region the fluxoids turn into 
a glassy flux phase (GFP). Hence, deep into the GFP the 
superconductor is still able to effectively shield external 
magnetic fields For the present sample, this can 

be deduced by the low-T saturation value x'ac — — l/47r 
clearly indicating that the superconducting phase is ex- 
tending over the total sample volume independently on 
i/dc- On the other hand, in the high-T region (still for 
T < Tc), mutual interactions and pinning mechanisms 
are no longer effective in comparison with thermal ener- 
gies so that the fluxoids are basically free to move, lead- 
ing to dissipation in turn. This is the so-called liquid flux 
phase (LFP). 

In correspondence with the sharp drop of Xac the imag- 
inary component Xac displays a single bell-shaped contri- 
bution peaked around a temperature value Tp also com- 
mon to the maximum in dXacl^T (see the inset of Fig. [1]). 
This results from the matching of t'm to the inverse cor- 
relation time 1 /tc characteristic of the motion of vortices 
through the crossover between GFP and LFP [13, [H [Ti- 
ll^. The most natural way to investigate in the vortex 
dynamics in detail is then to repeat Xac vs. T scans at 
different or, equivalently, to perform scans at dif- 
ferent T values across the GFP-LFP boundary. In the 
former case one extracts the Tp values from the position 
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FIG. 1: (Color online) Main panel: Xac vs. T data for both 
real and imaginary components at v^n ~ 9685 Hz and at dif- 
ferent Hdc (after correction for the demagnetization factor). 
Inset: zoom of data for Hdc = 30 kOe around the step-hke 
drop of the real component. The dashed line represents the 
derivative of xLc with respect to T (arb. units). 



of the maxima of Xac vs. T curves at different employ- 
ing empirical peak functions around the maxima (see the 
close symbols in the main panel of Fig. [2] and the rela- 
tive upper inset). In the latter case, one associates to the 
selected T the characteristic value where a maximum 
in Xac appears (open symbols in Fig. [2l raw data not 
shown). Both kinds of measurements were performed on 
the considered single crystal, always confirming a sizeable 
dependence of Tp on v^. Given the extreme narrowness 
of the observed transitions, anyway, even uncertainties in 
the thermalization up to ±0.01 K can be detrimental for 
the latter measurements, as indicated by the wide error 
bars in Fig. [2l Correspondingly, a comparison between 
T-swept curves was chosen as the most reliable method 
by always keeping under the same conditions the slow 
warming rate through the transition. 

If the degree of quenched disorder in the system is 
sizeable and, correspondingly, the density of pinning cen- 
ters is high, thermally-activated processes are the leading 
mechanisms in driving the dynamical relaxation of the 
fluxoids. This results in a logarithmic dependence of Tp 
on and is generally interpreted in terms of thermal^ 
activated flux creep (TAFC) of the vortices 0, Ei 



In the presence of extremely pure and clean supercon- 
ducting systems in the limit of weak disorder, on the 
other hand, the lower concentration of pinning centers 
makes the correlations and interaction among the vor- 
tices to clearly appear and to play a key role in the 
magnetic relaxation. This can allow, for instance, the 
detection of a scaling critical behaviour associated with 
a true thermodynamical phase transition between LFP 
and GFP. As it is well known, the fingerprint of the 
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FIG. 2: (Color online) Main panel: position of the maximum 
in Xac vs. T curves as a function of (full symbols). Open 
symbols refer to the estimates from isothermal Xac vs. 
curves. The dashed line shows the expected logarithmic trend 
typical of thermally-activated processes. Data are taken at 
-ffdc = 250 Oe (data for H^^ = 90 kOe are plotted in the 
lower inset with the same meaning of symbols). Upper inset: 
Xac vs. T at _ffdc = 250 Oe and at different v^- 



critical behaviour is the appearance of power-law func- 
tional forms describing the trend of physical quantities 
as a function of some reduced variable quantifjdng the 
distance of the system from the critical point [21[. Our 
results clearly reveal such a scaling behaviour instead of 
a logarithmic one in the measured single crystal. In the 
case of iJdc — 250 Oe shown in the main panel of Fig. [2l 
in fact, the departure from the dashed line (representing 
a guide to the eye for the expected trend emerging from 
a TAFC scenario) is extremely evident. The effect, even 
if less marked, is present also at higher values of the ap- 
plied static field, as shown in the lower inset of Fig. [2] in 
the case of iJdc = 90 kOe. Similarly to what is reported 
in the literature also in the case of cuprate superconduc- 
tors (see Ref. (3), a weak dependence of Tp on has 
been detected even at the nominal iJdc = Oe value (in 
the present sample, the variation of Tp at TJdc = Oe 
over the entire range of is less than 0.05 K). Neverthe- 
less, it is stressed that such effect should be considered 
as spurious mainly due to strong demagnetization effects 
(associated with the high value of Dm) possibly enhanc- 
ing the sum of both H^c and any residual component of 
Hdc to effective values close to Hd- 

With reference to the model discussed in Refs. [l|, [l^, 
one can assume that the transition of vortices from LFP 
to GFP occurs at a critical value Tg dependent on the 
value of the external static field. From dynamical scaling 
considerations the following relation can be derived for 
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FIG. 3: (Color online) Main panel: Tp vs. u-^n data at //dc = 
250 Oe (already plotted in the main panel of Fig. ^ after 
analysis according to Eq. (see continuous line). The data 
for Hdc ~ 90 kOe have been reported in the lower inset as in 
Fig. [2] (the error bars are within the symbols). Upper panel: 
dependence of the critical exponent on -ffdc- 

the reduced temperature eg(fm) 

eg(l^m) = 7f = — (1) 

where Tp^v^) is straightforwardly derived from the max- 
ima in Xac vs. T curves as discussed above. The param- 
eter Pq in this expression is a microscopic characteristic 
attempt frequency with a typical value in the THz range. 
The dimensionality of the system is indicated by D while 
z and v are the dynamical and static critical exponent 
for the vortex-glass correlation time and length, respec- 
tively. One can then deduce the values for Tg at the 
different values of H^c from plotting the data presented 
in the main panel of Fig. [2] according to Eq. ([Ij. The ac- 
curacy associated with the assessment of Tg values is re- 
markably high 0.01 K) since, by taking values beyond 
this range, results clearly deviate from the power-law be- 
haviours evidenced in Fig. [3] (both main panel and lower 
inset). This high precision in the estimate of Tg allows 
us to sizeably reduce the mutual statistical dependence 
among the fitting parameters in Eq. ^ . 

The results of the fitting procedure of the experimental 
data are in good agreement with theoretical expectations. 
The attempt frequency vq, in particular, is found to take 
values vq ^ 10^^ 10^^ Hz slightly higher than other re- 
ports in the literature but nonetheless reasonable. The 
dependence of the critical exponent on i?dc, on the other 
hand, is shown in the upper inset of Fig. [3] Most inter- 
estingly, a systematic increase can be discerned with in- 
creasing Hdc until the saturation value 7.8±0.2 is reached 
for Hdc ^ 20 kOe. This saturation value is in close agree- 
ment with previous reports on cuprate compounds like 
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YBasCusOr and BiaSraCaCuzOg+a [llii-lil where val- 
ues ~ 6.5 8 are typically reported. On the other hand, 
the large change evidenced in the low-field regime cannot 
be explained in terms of a crossover between different val- 
ues of dimensionality for the vortex phase. Following the 
arguments reported in Ref. |24| . in fact, one can estimate 
the typical magnetic field inducing such dimensionality 
crossover as Her ~ 4<i>o/(i^7^. For the inter-plane dis- 
tance one can assume d ~ 6.5 A (see, for instance, Ref. 
and references therein) while 7 ~ 2 holds '25!|, result- 
ing in Hci values as high as 5 x 10^ Oe. This is similar 
to the case of 1111 compounds even if, in that case, the 
field dependence of the critical exponent is the opposite 
of what is observed here Q . Our results instead closely 
resemble what is reported in the case of twinned crystals 
of YBa2Cu307_5 where the field-dependence of the criti- 
cal exponent typically shows sudden step-like increase at 
some characteristic field values of the order of few Teslas 
[iil [53] ■ The interpretation of the results in those cases 
are typically given in terms of a crossover between dif- 
ferent types of glass phases in relation with the typical 
dimensionality of crystalline defects. Further studies on 
the present single crystals are needed in order to gain 
more insights into the detailed features of vortex pin- 
ning and their connection to the resulting glassy phase. 
In this respect, as discussed in previous reports about 
MgB2 samples [28|, it should be anyway stressed that 
low values of magnetic field [Hdc ^ 10 kOe) are expected 
to favour the emergence of a Bragg-like glassy phase in 
presence of weak degrees of quenched disorder [29j and, 
in turn, to lead to a decrease of the values of critical ex- 
ponents. This well agrees with our results and strongly 
supports the same scenario in Ba(Feo.9Coo.i)2As2. The 
high-field region characterized hj v {z + 2 — D) ^ 7.8 is 
then associated to a glass state with a higher degree of 
amorphousness (see Fig. |4|). 

The Tg values estimated from the fitting procedure ac- 
cording to Eq. ([T]) are plotted as a function of the mag- 
netic field in Fig. 21 As reported in the same graph, their 
dependence on the corresponding H^c — Hg values can 
be nicely fit by a mean- field like function [rc(0) — Tg] oc 
(i?g)^^^ with the fitting parameter /3 — 1.30 ± 0.05 in 
excellent agreement with the expected result /? = 4/3 
from the vortex-glass model for Z) = 3 and with previous 
reports in YBa2Cua07 [H, [111- According to the model 
discussed in Ref. |29j, further measurements of magnetic 
ac susceptibility at higher values of magnetic field allow- 
ing to detect a critical endpoint of this transition line 
could possibly further check the validity of the Bragg- 
glass scenario described above. Results for Tg (iJdc) are 
completed by the T-dependence of H^2 in Fig. ID Ex- 
perimental points have been derived from measurements 
of FC magnetization at different H^c values by means 
of a Quantum Design vibrating-sample SQUID magne- 
tometer (the fitting procedure of the diamagnetic onset 
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FIG. 4: (Color online) Phase diagram for the superconducting 
state of Ba(Feo.9Coo.i)2As2. The continuous line is a fit to 
the Hg data according to the mean-field function discussed in 
the text. The shaded areas refer to the two different regimes 
for the dependence of the critical exponent on H^c presented 
in the upper inset of Fig. [3] and discussed in the text. The 
dashed line is a linear fit to the data (see text). 



has been described previously in Ref. 11). The relative 
raw data are not shown. A linear fitting of the data 
leads to a slope value —2.50 ±0.05 T/K, in perfect agree- 
ment with previous reports on different single crystals of 
Ba(Fei_^Co:r)2As2 

In summary, the extremely high quality of the 
investigated single crystal of superconducting 
Ba(Feo.9Coo.i)2As2 allows us to evidence by means 
of magnetic ac susceptibility measurements a critical 
behaviour associated with the phase transition of vortex 
lines from liquid to glassy states. Results in the high- 
field regime are in quantitative agreement with typical 
reports for cuprate superconductors, suggesting a shared 
underlying physics for both kinds of high-Tc supercon- 
ductors. The field-dependence of the critical exponent 
cannot be explained in terms of any dimensionality 
crossover for the vortex phase. Instead, the results 
suggest a non-trivial transition among different types 
of glassy-states, where the degree of amorphousness 
is decreased towards the limit of Bragg-glass with 
decreasing the magnetic field. Our results open the 
way to a more detailed investigation of the physics of 
vortices across the phase diagram of Ba(Fei_a;Co2;)2As2. 
The evolution of the degree of glassiness and its possible 
interplay with magnetic fluctuations arising from the 
progressively enhanced nesting of the Fermi surface with 
reducing x, in particular, could be an extremely peculiar 
study based on our present data. 
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